Andreev Conductance of Chaotic and Integrable Quantum Dots 
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We examine the voltage V and magnetic field B dependent Andreev conductance of a chaotic 
quantum dot coupled via point contacts to a normal metal and a superconductor. In the case where 
the contact to the superconductor dominates, we find that the conductance is consistent with the dot 
itself behaving as a superconductor- it appears as though Andreev reflections are occurring locally 
at the interface between the normal lead and the dot. This is contrasted against the behaviour 
of an integrable dot, where for a similar strong coupling to the superconductor, no such effect is 
seen. The voltage dependence of the Andreev conductance thus provides an extremely pronounced 
quantum signature of the nature of the dot's classical dynamics. For the chaotic dot, we also study 
non-monotonic re-entrance effects which occur in both V and B. 



X3 

o 



> 

o 
o 
o 



I 

o 
o 



X 



I. INTRODUCTION 

Though the actual process of Andreev reflection is sim- 
ple to describe - an electron in a normal metal incident 
on a superconductor is reflected back as a hole -, 
it serves as the fundamental basis for some of the most 
striking effects knovifn in mesoscopic physics In par- 
ticular, Andreev reflection may be viewed as the phe- 
nomenon underlying the proximity effect, in which a su- 
perconductor is able to strongly influence the properties 
of a nearby normal metal. Andreev reflection is also the 
process responsible for the unique conductance proper- 
ties of normal-superconducting (NS) junctions, such as 
the conductance enhancement observed in NS point con- 
tacts. 

Attention has recently turned to so-called "Andreev 
billiards" as ideal systems in which to study the prox- 
imity effect P-p^. Such systems consist of an isolated 
normal metal region, small enough that electrons remain 
phase coherent within it (i.e. a quantum dot), coupled 
weakly through a point contact to a superconducting 
electrode. Calculations of the density of states in such 
structures have shown, rather remarkably, that the prox- 
imity effect is sensitive to the nature of the classical 
dynamics- a gap of size A in the spectrum is present 
in the case of a chaotic billiard, whereas in the integrable 
case the density of states vanishes linearly at the Fermi 
energy [|§. 

In the present work we also focus on the Andreev bil- 
liard, but now add a second point contact leading to a 
normal electrode and investigate the Andreev conduc- 
tance of the resulting structure. The Andreev conduc- 
tance is the differential conductance dl/dV at voltages 
smaller than the superconducting gap in the S electrode, 
where Andreev reflection at the interface between the 
normal metal dot and the superconductor is the only cur- 
rent carrying mechanism. Unlike previous studies |]l9| , ^ , 
we calculate the full voltage {V) and magnetic field {B) 
dependence of the conductance. 

Two questions are of particular interest in this study. 



First, does the sensitivity of the proximity effect to 
chaotic versus integrable dynamics, as seen in the den- 
sity of states, also manifest itself in the Andreev conduc- 
tance? We find that indeed it does. In the case where the 
contact to the superconducting electrode is much wider 
than that to the normal electrode, the V and B depen- 
dent conductance of a chaotic dot is consistent with as- 
suming the dot has itself become superconducting. The 
conductance is doubled with respect to the normal state, 
and remains^ bias-independent until eV reaches the in- 
duced gap A in the Andreev billiard. This is in sharp 
contrast to what is found in exact calculations for inte- 
grable dots. For both rectangular and circular dot ge- 
ometries, the conductance drops off linearly with voltage 
without any plateau near V = 0. We argue that this 
difference is ubiquitous for integrable versus chaotic sys- 
tems. 

The second motivation for the present study is the 
question of re-entrance effects- is the behaviour of the 
conductance simply monotonic in V and B7 In the case 
of a diffusive NS junction, it is well known that this is 
not the case [ ]lT| , p^ . At zero voltage and also at volt- 
ages large enough to break electron-hole degeneracy, the 
conductance of the junction is the same as in the normal 
state. However, a conductance enhancement does occur 
at intermediate voltages. We examine re-entrance effects 
in the conductance of chaotic Andreev billiards, both in 
the case of ballistic contacts, and in the case where both 
point contacts contain opaque tunnel barriers. 

The remainder of this paper is organized as follows. In 
Section II, we present our model for the chaotic dot and 
introduce the technique used to calculate the voltage and 
field-dependent Andreev conductance. In Section III, we 
discuss results in the case where the contact to the super- 
conductor dominates the contact to the normal electrode, 
and contrast the results for a chaotic billiard with exact 
quantum mechanical calculations for two integrable sys- 
tems. In Section IV, we discuss re-entrance effects. We 
conclude with a synopsis of our key results in Section V. 
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II. MODEL AND TECHNICAL DETAILS 



A. Formulation of the Problem 




FIG. |l[ Schematic drawing of a dot coupled via point con- 
tacts to a normal metal and a superconductor. 

We consider a chaotic quantum dot coupled via point 
contacts to a normal metal and a superconductor, each 
having respectively and Ns propagating modes at 
the Fermi energy Ep (see Fig. |l|). We assume that the er- 
godic time is much shorter than other relevant timescales 
of the dot (i.e. the inverse superconducting gap h/A and 
the dwell time), so that random matrix theory (RMT) 
p3[ may be used to describe its transport and spectral 
properties. In RMT, the Hamiltonian of dot is repre- 
sented by a M X M Hermitian matrix H which, at zero 
magnetic field, is real symmetric and a member of the 
Gaussian orthogonal ensemble: 



P{H) = exp(--A/A"^tri/2). 



(2.1) 



The matrix size M is sent to infinity at the end of the cal- 
culation. The energy scale A is related to the mean level 
spacing 25 by A = More specifically, 25 is the mean 

level spacing for particle-like excitations in the absence 
of a coupling to the superconductor; with the supercon- 
ductor, the relevant excitations are of mixed particle-holc 
type, and have a level spacing 6 for excitation energies 
much larger than the gap energy. 

In the case of a non-zero magnetic field, the Hermitian 
matrix if is a member of the Pandey-Mehta distribution 
00: 



(2.2) 



As 7 increases from to 1, the distribution evolves from 
one with complete time-reversal symmetry (Gaussian or- 
thogonal ensemble) to one with no time-reversal sym- 
metry (Gaussian unitary ensemble). The parameter 7 
may be related to the magnetic flux <I> through a two- 
dimensional dot having area A |lq,|6[ : 



M-f^ = C 



$0 



AS 



(2.3) 



where $0 = hc/e is the flux quantum, vp is the Fermi 
velocity, and C is a constant of order unity. 



The normal (TVjv + Ns) x {Nn + Ns) scattering ma- 
trix S{s) of the system at an energy e above Ep can be 
expressed in terms of the matrix H [ [TtI] : 

S{e) ^ 1 - 2mW'^ {e - H - inVVW"')'^ W (2.4) 

rNNie) tNsis) 
tswie) rss{e) 

where W is an M x {N^ + Ns) matrix representing the 
coupling between the point contacts and the dot, having 
elements: 



W„ 



6^n-i2MS) 



1/2 
1/2 



2 - y„ - 2^1^ 



1/2 



(2.5) 



The r„ are the transmission probabilities for each mode, 
which we take simply as Tn for modes coupled to the 
normal electrode, and Ts for those coupled to the super- 
conducting electrode (this in turn defines Zn and Zs). 

For voltages below the excitation gap A of the su- 
perconductor, electrons and holes incident on the dot- 
superconductor interface may be Andreev reflected. In 
this case scattering from the dot, as seen from the nor- 
mal metal contact, can be represented by a 2Nn x 2N]\[ 
scattering matrix S: 



S = 



ree{e) reh{e) 
rheis) rhhis) 



(2.6) 



Here, r^eis) is the N^ x Nn matrix describing the An- 
dreev reflection of an incoming electron in the N point 
contact to an outgoing hole in the same lead; rehis), 
ree(e) and rhh{s) are defined analogously. These ma- 
trices may be written in terms of the sub-matrices of the 
normal scattering matrices S{e) B. For Vheie) one has: 



with 



rhe{e) = i^s(-e)Mee(e)a(e)isAr(£), (2.7) 

M,,{s) = [{1 - aie)\ssie)r*ssi-e)r' 
a{e) — exp(— i arccos(£/A)) 



For voltages below the gap A in the S electrode, the 
zero temperature conductance of the system is given by 
the Tabikane-Ebisawa formula 11131 : 



G(eV) = ^tr rl{eV)r^^ieV). 



(2.8) 



We wish to calculate the ensemble-averaged Andreev 
conductance of the dot for arbitrary values of voltage 
V and magnetic field B; previous studies focused exclu- 
sively on the cases where V, B were or large enough to 
completely break the symmetry between electrons and 
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holes. To this end, we first rewrite Eq. ( p.6[ ) in a manner 
which is formally equivalent to the normal-state expres- 
sion (U). Letting Qs ^ nWsWl and f^^v = ttWnWI, 



we define the 2M x 2M effective particle-hole Hamilto- 
nian 



n = 



H 






H* 



(2.9) 



and the self-energies from the leads: 




A 







VA2 - £2 



(e/A)05 ns 
ns (e/A)Os 



(2.10) 



With these definitions in place, a direct algebraic manip- 
ulation shows that the particle-hole scattering matrix S 
in can be expressed in terms of the effective retarded 
Green function Qn, 



by 

S{e) = \-2m 



K 

w\ 



N 



Wn 
Wn 



(2.11) 



(2.12) 



The major simplification offered by Eq. (2.12) com- 
pared to Eq. (2.6) is that it allows one to compute the 
conductance in a manner analogous to that used in the 
normal state, as we now demonstrate. From this point 
onward, we focus on the case A ^ e, eV . In this regime, 
the properties of the system become independent of the 
specific details of the superconductor. Calculations re- 
taining a finite A will be presented elsewhere ||2Ct|. 



B. Details of the Calculation - B = Q 



We proceed to average Eq. ( p.8D for the Andreev con- 
ductance over the Gaussian orthogonal ensemble defined 



by Eq. (2.1), which is the appropriate ensemble for zero 
magnetic field. We make use of the relation 



\2 

(HijHki) = {SikSji + 6ii6jk)—- 



(2.13) 



By expressing the scattering matrix in terms of the Green 
function Q, the trace appearing in the conductance for- 
mula may be represented as a standard "conductance 
bubble" diagram (see Fig. Further, at this stage 

of the calculation }? /M is taken to be a small parame- 
ter, meaning that the usual diagrammatic technique for 
impurity averaging may be used. To obtain the leading- 
order result in l/7Vjv,l/A^s, we need only sum diagrams 
which have no crossed lines (lines indicate the averaging 
of matrix elements of if using the rule ( 2 . 1 E ) ; see Fig . ||) . 



Note that A will be sent to infinity only at the very end of 
the calculation; this procedure corresponds to first tak- 
ing the matrix size to infinity while keeping a finite band- 
width, and then independently taking the bandwidth to 
infinity. This ordering of limits is necessary to generate 
a well-defined perturbative expansion. 



(m,i)*(n,i) 



(m-,i-) ^ {n-X) 



L J(m,i) (n.i'j 



d) (m.i) 



=i<(m,i-) 



FIG. ^. (a) Definition and diagrammatic representation 
of the unaveraged Green function. This is a 2M x 2M ma- 
trix; n,m refer to mode indices, while are electron-hole 
indices which can take one of two values (e or h). (b) The 
effective Hamiltonian 7i is represented by a cross; the dashed 
lines indicates averaging, (c) Diagrammatic Dyson equation; 
the double line is the averaged Green function (5). (d) Self 
energy to leading order in 1/M; the neglected higher order 
terms correspond to graphs with crossed dashed lines. The 
intermediate index m is to be summed over. Note that E and 
hence (£/) are diagonal in mode space but have oflF-diagonal 
terms in particle-hole space. 
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FIG. (a) direct contribution to the Andreev conduc- 
tance. The trace appearing in Eq. (2.8) is represented by the 



bubble diagram. The double lines indicate the averaged Green 
function {Q). The index n corresponds to a mode coupled to 
the normal contact and should be summed over, (b) Dif- 
fusion contribution to the Andreev conductance, (c) Dyson 
equation for the AM x 4M matrix D. Note that upper and 
lower branches need to have matching mode indices, but not 
matching particle-hole indices. 
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The first step in this framework is to calculate the av- 
eraged matrix Green function {Q). As in 1^,^, we find 
the following self-consistent Dyson equation: 

{g{e)) = [e- Y.l{e) - (e) - ' , (2.14a) 

where the self-energy from averaging E is given by 

Eee Ee/i 



S = 1m «> 



1 



M ' 



MK-itrGhe) (trGhh) 



(2.14b) 



In this equation, 1m denotes the Af x M unit matrix, 
and iS> denotes a direct product; it is pictured diagram- 
matically in Fig. (^). Having performed the summation 
of diagrams, we now let M ^ oo keeping £ , (5 , Ns and 
iVjv fixed. 

Three of the equations thus obtained relate the com- 
ponents of the self energy to one another: 



-'ee — ^hh j ^eh — ^he: 



E2 \ 2 



(2.15a) 



(2.15b) 



The last equation allows us to parameterize S in terms 
of a pairing angle 0{e): 



Asin6'(e) = -T,eh{£), 
\cos0{e) = il]ee(e)- 



(2.16a) 
(2.16b) 



The remaining self-energy equations now take the form: 



tan(6'(£)) 
Qn 



NsQs 



^(2-rA.[l-COS0(£)])-\ 



)s='^{2-Ts[l + sme{e)])-\ 



(2.17a) 

(2.17b) 
(2.17c) 



To obtain a unique solution of the self-energy equations 
we have imposed the boundary condition See(e) ^ —iX 
as £ ^ CXI ; this represents the physical condition that we 
recover a normal metal with a constant density of states 
at large energies. The pairing angle 9 can is related the 
density of states of the dot: 

p(£)^--tr(g(£ + 0+)) = ^^[^"f("»] ; (2.18) 

9 ^ 71/2 corresponds to a fully superconducting state, 
while 9^0 corresponds to the normal state. It is inter- 
esting to note the analogy between 9(e) and the pairing 



angle used in the quasi-classical theory of dirty normal- 
superconducting interface s [p2[ and in the circuit theory 
of Andreev conductance |23|] . In the present case, the 



normalization condition (2.15b) does not need to be ex- 
ternally imposed, but is a direct consequence of the av- 
eraging procedure. 

While the leading order solution of the self energy S is 
sufficient if one is only interested in the density of states 
(as in ref. [^), the conductance calculation requires that 
Ns/M , Nn/M and e/MS corrections to Eq. ( ^.15bD be 
calculated. This correction may be expressed in terms of 
the leading order self energy solution: 

- 1 - - (nsQs (sin(g) + Zs) (2.19) 



+ NnQn (cos(6I) + Zn) + i 



Having computed the self-energy from averaging S 
and thus the average Green function (G), we can pro- 
ceed to sum diagrams for the conductance. In analogy 
to the usual impurity technique [^4| , two sets of con- 
tributions arise (see Fig. ^a) to leading order: G = 
(4e^//i)((7Dir + ffDiff)- The first is a direct contribution 
which is completely determined by (G), that is, by the 
ensemble- averaged probability amplitude {rhe)- 



ffDir = (2^)'tr W],{G)W^WI{G)^W^ 



' N 



= tr {rhe){rhey 



(2.20) 



This contribution may be interpreted as arising from An- 
dreev reflections which effectively occur at the interface 
of the normal lead and the cavity. The second contribu- 
tion is due to the fluctuations of r^^: 



5Diff 



{reh - (reh)) {reh ~ {reh)y 



(2.21) 



It describes the current carried by quasiparticles in the 
dot which are Andreev reflected at the interface of the 
dot and the superconductor. Diagrammatically, it is 
equivalent to a diffusion ladder, where the averaging 
links the upper and lower branches of the conductance 
bubble (see Fig. |b-d) We ignore here quantum 

corrections, which are formally smaller by a factor of 
max{l/N]y,l/Ng). Note that in computing the diffu- 
sion sum, each individual graph is of order 1/Af; this 
means that terms of order 1/Af must be retained when 
computing the 4x4 matrix inverse arising from summing 
the series. As with the calculation of the average Green 
function (G), we only let A/ tend to infinity after having 
performed the partial summation of diagrams. 

To present the results of the conductance calculation, 
we first define the following kernel functions: 



J 



4 





Qn 







|l + ZArC0s(6')|^ +27Vs 1 



Qs 



|l + Zssin(0)r 



2Rg 



NnQn {Zn + cos(0)) + NsQs {Zs + sin(0)) + — Im [cos(0)], 





7re, 



= 2Ar^,|g^,sin(6')|V2Ars|(3scos(6l)|^ 



nw(£) = 2A^Ar|Qjvp{ (1 + 6^2, + Z%) (1 + I cos((?)|2) + SZat (1 + Zl) Re [cos(0)]}, 
ns(e) = 2Ns\Qs?[{l + ^Zl + Z|) (l + | cos(f?)|2) + 8Z5 (l + Z|) Re [cos(0)]}, 



(2.22a) 

(2.22b) 
(2.22c) 

(2.22d) 



r 



where Zjq and Z5 were defined below Eq. (2.5). With 
these definitions, we find: 



2iV, 



N- 



rpl 

^N 



|2-2TArsin2(6l/2)| 



|sin(0)r 



(2.23a) 



5Diff = 8A^; 



(a2 - r!^ 



nAr|sin(6')p (2.23b) 



(l-Z2,)'A|sin(0)p+ns|cos(0)P 



Equations ( p.22aD - ( p.23b| ), along with Eq.(pl^) for the 
self energy S, determine the ensemble-averaged i? = 
Andreev conductance through the dot to leading order 
in l/A^AT , l/7Vs for all voltages such that < A. The 
extension of these formulas to non-zero B are presented 
in the Appendix. 



III. PROBING INDUCED 
SUPERCONDUCTIVITY VIA CONDUCTANCE 

A. Chaotic Dot 

In this section, we consider the situation where the cou- 
pling between the dot and the superconductor is much 
stronger than the coupling between the dot and the nor- 
mal lead, NnTn ^ NsTs- In this limit, the normal 
metal will only weakly perturb the properties of the dot- 
superconductor system; thus, we may view the conduc- 
tance through the structure as being a probe of the in- 
duced superconductivity in the dot. As mentioned in the 
Introduction, previous studies have examined the den- 
sity of states of such Andreev billiards in the absence of 
a normal lead |^ . It was found for the case of a chaotic 
dot that an energy gap A opens up on a scale set by the 
inverse of the time needed for a particle to escape from 
the dot to the superconductor: 



A = c{Es) 
, h 

Es = 



TS., 



NsTsS 
2tt 



(3.1) 
(3.2) 



The parameter c is of order unity and a monotonic func- 
tion of T5; it varies from 0.6 in the case of no tunnel 
barrier, to 1 in the case of an opaque tunnel barrier. It 
was also found that the shape of the density of states 
above A was vastly different in these two limits: in the 
tunnel regime ^ 1, the density of states was BCS- 
like, having a square-root singularity at A, whereas in 
the case of no tunnel barrier, the density of states grad- 
ually increased above A. Note that the gap A does not 
depend on the gap A in the bulk superconductor, which 
was taken to infinity in the calculations. 

The question that naturally arises here is how the 
induced superconductivity seen in the density of states 
manifests itself in the conductance. Some insig ht may be 
obtained by considering equations ( ^.23a ) and ( 2.23b ) to 
lowest order in the small parameter 77 = N^Tn/NsTs- 
Without the normal lead, we have a gap in the dot den- 
sity of states up to an energy A. As the density of 
states is proportional to Re [008(0 (e))], this implies that 
Re [e{e)] = 7r/2 -I- 0(7/) for e < A . Using this, the "sub- 
gap" (i.e. £ < A) direct contribution to the conductance 
takes the form: 



9 = 5Dir(e) 



2NmT] 



nJ-n 



(2 - T^)2 - 4(1 - rw)(Im[sin(0o(e))])2 
+ 0{7j), (3.3) 



while the diffusion contribution goiff is of order rj^ and 
thus negligible. As the energy is increased above the gap, 
the density of states returns to its normal-state value, and 
consequently 9{e) ^ to leading order. For e A the 
direct term is negligible (being proportional to sin(0)), 
while the diffusion term gives the normal-state conduc- 
tance: 



5Diff( 



NmT, 



N-^N 



oiv). 



(3.4) 



The above considerations become extremely suggestive 
when one considers the case T/v = 1- There is a perfect 
conductance doubling for voltages below the effective gap 
A, while at higher voltages the conductance drops to its 
normal state value: 



2iVAr, 

Nn, 



for e < A 
for e > A. 



(3.5) 
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FIG. ^. Conductance versus voltage for Ns = IO^Nat, 
Tat = 0.1, illustrating the "local Andreev reflection" effect, 
in which the dot itself acts as a superconductor having a gap 
A. Solid lines are calculated curves for, from top to bottom 
Tn = 1, 0.9, 0.75. The dashed curve shows the BTK resuh for 
Tm = 1.0; close agreement is found for all values of Tn- The 
dot density of states (shown in inset) in this case resembles 
a BCS density of states. The voltage is measured in units of 
Es, the inverse of the escape time to the superconductor (cf. 
Eq. Q). 
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FIG. Q. Conductance versus voltage for A''s = 10 Nn, 
Ts — 1, illustrating the "local Andreev reflection" effect. 
Solid lines are calculated curves for, from top to bottom 
Tn ~ 1,0.9,0.75; dashed curves show BTK result for same 
values of Tat. The differences in the dot density of states 
(shown in inset) from a BCS density of states leads to large 
differences from the BTK lineshapes. 

This is precisely what would be expected if the nor- 
mal point contact were in perfect contact with a bulk 
superconductor having a gap A- conductance doubling 



would be expected below the gap due to Andreev reflec- 
tion, whereas above the gap the normal-state conduc- 
tance would be recovered as quasiparticles would now 
carry the current. In this respect, note that below the 
induced gap A both the averaged reflection probability 
(kehP) and the averaged reflection ampHtude |(reft)| are 
non-zero; for eV ^ A, the latter vanishes. Nonetheless, 
in both cases (i.e. for voltages above and below A) An- 
dreev reflection is the only current carrying process. For 
energies below A, the dot itself appears as if it were su- 
perconducting and Andreev reflection effectively occurs 
at the normal lead-dot interface. For energies far above 
A, current through the dot is effectively carried by quasi- 
particles in the dot (which are either electron or hole-like) 
which Andreev reflect at the dot-superconductor inter- 
face. 

The quantum dot continues to act as if it were a su- 
perconductor in the case where the contact to the nor- 
mal lead is no longer perfect (i.e. ^ 1); one need 
only note that if Im[sin(6'o(e))] is replaced by e/A, Equa- 
tion ( |3.3| ) becomes identical to the Blonder-Tinkham- 
Klapwijk (BTK) formula for the sub-gap conductance of 
a normal point contact - superconductor junction, where 
there is a tunnel barrier at the interface having transmis- 
sion Tn ■ The BTK formula is derived assuming there 
is no spatial separation between the sites of Andreev scat- 
tering and normal scattering; the only energy dependence 
thus comes from the Andreev reflection phase a{e), not 
from the lack of electron hole degeneracy at finite volt- 
ages. The replacement of e/A by Im[sin(^?o(e))] in Eq. 
( |3.3|) means that the Andreev reflection phase a in the 
present case is not necessarily the usual a{e) given in Eq. 
( ^.7|) , but is rather defined through: 

1 + Re[5(£)2] 



lm[sin(0o(£))]' = 



(3.6) 



This notwithstanding, the overall implication is still that 
for eV < A, we can effectively consider all Andreev re- 
flections as occurring at the N-dot interface, as opposed 
to at the dot-S interface- the dot is indeed acting as 

though it were itself a superconductor. 

We have numerically solved Equations ( 2.17a| )- (|2.17c|) 
for the pairing angle 9(e), and used this to compute the 
Andreev conductance versus voltage. In Figures |and|, 
we plot the calculated conductance vs. voltage curves for 
various values of Tn , and compare to what would be ex- 
pected from the BTK theory for a simple N-S interface. 
We find an excellent agreement in the case of T5 ^ 1 
(see Fig. ^ even for voltages above the gap A in the 
dot. Such an agreement is not unlikely, as for Ts <C 1 
the dot density of states is BCS-like, and consequently 
the effective Andreev phase a is just equal to the usual 
Andreev phase a. In the opposite case of a transparent 
contact to the superconductor (Fig. clear deviations 
from the BTK lineshapes are seen. These deviations re- 
sult completely from the fact a ^ a, which is to be ex- 
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pected as the dot density of states in this case is quite 
different from the BCS form. 

The "induced superconductivity" effect in the dot also 
manifests itself in other manners. A straightforward 
calculation in the limit Tm shows that the An- 
dreev conductance becomes proportional to the dot den- 
sity of states (i.e.Re[cos(6')]); thus, the Andreev conduc- 
tance becomes equivalent to a conventional superconduc- 
tor tunneling-density of states measurement. We have 
also calculated the magnetic field dependence of the An- 
dreev conductance in the limit of small rj (see Appendix 
A); here too, our results are consistent with a picture in 
which the dot itself acts as a superconductor. The con- 
ductance enhancement at field remains constant until 
a critical flux given by: 



(3.7) 



where r^rg is the ergodic time, and C is a geometry- 
dependent constant of order unity. This is the same crit- 
ical field required to close the gap A. Note that unlike 
a conventional BCS superconductor, where the critical 
field is proportional to the gap, in the present case, Eq. 

A. 



[rhejs.c. 



dLe~P{L) = 



1 - i 




(3.7) implies that the critical field is proportional to 

The conclusions reached here are markedly different 
from what would be expected from a naive trajectory- 
based semi-classical analysis Consider the situa- 
tion of a ballistic dot with no tunnel barriers, where 
Ns » Nfq . In the semiclassical picture, electrons en- 
tering the dot from the normal lead typically bounce off 
the walls of the dot several times before hitting the su- 
perconductor, where they Andreev reflect. Ai V — 
holes are the time-reversed partners of electrons; thus 
when an Andreev reflection occurs, the hole will retrace 
the path of the incoming electron and cancel its acquired 
phase. Andreev trajectories will thus interfere construc- 
tively, leading to a large conductance enhancement at 
V ~ Q and a non-vanishing average reflection amplitude 
(rfie)- This enhancement should be lost however as the 
voltage is increased- at finite V, electrons and holes are 
no longer degenerate, and the phase acquired by the hole 
will not precisely cancel that acquired by the electron. 
The presence of a residual phase S(j) = eVL/ (Tivp) (where 
L is the length of the trajectory) will lead to destructive 
interference, and a consequent decrease in the conduc- 
tance. In this picture, even a small voltage should impair 
the conductance enhancement seen at zero voltage. 

To make this picture quantitative, we may use the 
fact that our exact random matrix theory indicates that 
the conductance below A is proportional to \{rhe)\'^ (see 



Eq. (3.3)). Using the semiclassical Andreev phase 
S(f) — eVL/{hvp) and the fact that path lengths in 
chaotic systems have an exponential distribution func- 
tion P{L) = exjp{—L/L)/L we estimate the average 
of the semiclassical reflection amplitude as: 



eVL 



(3.8) 



where L is the mean path length to the superconductor; 
for a ballistic dot, we have L = vpTs.esc- Hence: 



(3.9) 



where we have used A = ch/rs.esc. (see Eq. (pT|)) with 
c ~ 0.6 in the case of a perfect contact to the supercon- 
ductor. The semiclassical approach thus predicts that 
|(^/ie)|^ (and hence the Andreev conductance) will fall 
off with voltage as a Lorentzian for eV < A. This is 
clearly at odds with our fully quantum mechanic calcu- 
lation at Tn = 1, which shows that |(r;ie)|^ — 1 for all 
voltages below A. We return to the discrepancy between 
the semiclassical and quantum mechanical calculations in 
the next subsection. 



B. Integrable Dot 




eV/E„ 



FIG. ^ Conductance (scaled by Nn) vs. voltage for a 
rectangular dot and a circular dot, each with Nm = 10, 
Ns ~ 30Nn (shapes shown in inset, with the size of the nor- 
mal lead exaggerated). Each plot was produced by averaging 
over small variations in Ep. Unlike a chaotic dot (dashed 
line), the Andreev conductance of the integrable dots has no 
"plateau" for eV < A. The conductance of the rectangle 
drops below due to a weak localization correction. The 



energy scale Es is defined in Eq. (3.2) 



In this subsection, we examine the conductance of an 
N-dot-S system where the classical dynamics of the dot 
are integrable. Previous studies [0,0 have indicated that 
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the proximity effect in the density of states is very differ- 
ent for chaotic and integrable biUiards; in the latter, there 
is no induced gap, but rather the density of states tends 
to zero hnearly at Ep. We argue here that the proxim- 
ity effect's strong sensitivity to chaos versus integrabiUty 
also manifests itself in the Andreev conductance. 

We consider two different integrable systems: a rect- 
angular dot and a circular dot. In each case, the dot 
is coupled via narrow leads to both a normal metal and 
a superconductor; there are no tunnel barriers. As in 
the previous subsection, we consider the case where the 
width of the normal contact is much smaller than the 
width of the superconducting contact (see inset of Fig. 

so that the normal contact serves as a probe of the 
proximity effect in the quantum dot. The scattering ma- 
trix of the system is computed by numerically matching 
wavefunctions across the structure, and the conductance 
then follows from equation ( |2.8| ). The number of modes 
in the normal lead was fixed at A^v = int {tpW/ir) [p7| ; 
we averaged results over small variations in kp which did 
not change N^. Our results are displayed in Fig. ^; sim- 
ilar results are obtained if one changes the position of the 
two point contacts. 

The results for the two integrable systems are simi- 
lar to one another, and differ significantly from what was 
found for a chaotic dot- as opposed to flat sub-gap region 
followed by a rapid drop-off, we have instead a gradual 
linear reduction of the conductance with voltage. We 
argue that the difference exhibited here is generic for in- 
tegrable systems. The reason is the same as that given in 
to explain the difference in density of states, namely, 
that the distribution of path lengths is very different for 
integrable and chaotic dynamics In the first case, 

there is a power law distribution of path lengths, mean- 
ing that there is an appreciable probability to find long 
paths. Even for small voltages, these long paths will 
quickly acquire a random phase, leading to destructive 
interference and a reduction of the conductance. In con- 
trast, chaotic systems have an exponential distribution 
of path lengths- long paths are quite rare. In this case, 
small voltages will not be able to cause any significant 
phase-randomization, and thus there will be no resulting 
destructive interference of Andreev trajectories. 

Though the semiclassical reasoning in terms of path 
lengths presented here provides a qualitative account for 
the difference between the chaotic and integrable An- 
dreev conductance lineshapes, attempts to translate it 
into a quantitative theory have not been successful. As 
was demonstrated in the previous section, a simple semi- 
classical theory for the chaotic case fails to recover the 
correct lineshape. A similar problem is encountered when 
one tries to do the analogous calculation for the inte- 
grable case- such a calculation predicts the Andreev con- 
ductance should fall off quadratically with voltage, albeit 
at a faster rate than in the chaotic case. This behaviour 
is clearly at odds with our exact calculation showing a 



linear dependence on voltage. In many ways, the failure 
of a semiclassical approach is not surprising. It is known 
that semiclassical approximations are unreliable in super- 
conducting systems, as the usual diagonal approximation 
is worse than for normal systems . 

Finally, we note that there is a striking connection be- 
tween the Andreev conductance's sensitivity to the na- 
ture of the dot's classical dynamics and weak localization. 
It has been shown both theoretically and experimen- 
tally that the magnetic field dependence of the weak 
localization correction of a quantum dot with two normal 
leads is very different for chaotic versus integrable dots. 
In the chaotic case, one finds a smooth Lorentzian field 
dependence, while in the integrable case a much sharper 
profile is found, with a cusp at zero magnetic field. The 
similarity between the Andreev conductance effect and 
that in weak localization is not coincidental; both effects 
rely on the interference of time reversed paths. 

Despite this strong similarity, it is worth noting that 
the effect in the Andreev conductance is much more pro- 
nounced. Here, the difference between the chaotic and 
integrable lineshapes is more severe than in the weak lo- 
calization effect (i.e. the chaotic case is much flatter than 
a Lorentzian). The magnitude of the effect is also much 
larger in the Andreev case; the signature of chaos versus 
integrability is the entire conductance itself, not a quan- 
tum correction like weak localization. For this reason, 
the effect in the Andreev conductance should be observ- 
able in a single sample, whereas an ensemble average is 
required in the weak localization case. 

IV. RE-ENTRANCE EFFECTS 

We shift focus in this section, and examine so-called 
"re-entrance" phenomena in the Andreev conductance of 
a chaotic quantum dot. These effects are loosely deflned 
by non-monotonic behaviour of the conductance in ei- 
ther voltage or magnetic field. They are well known in 
the case of diffusive NS systems , where the word 

"re-entrant" is used because the Andreev conductance 
is the same as the normal conductance at zero voltage 
and magnetic field, and at high voltage or field, but not 
for intermediate values. The theory developed here al- 
lows us to address this behaviour using a scattering ap- 
proach, whereas previous approaches almost exclusively 
made use of the quasi-classical Green function technique. 
In what follows, we discuss the cases of ballistic contacts 
(T/v — Ts = 1) and tunnel contacts {Tn,Ts <C 1) sepa- 
rately. 

A. Ballistic Contacts 

In the absence of tunnel junctions, the equations deter- 
mining the conductance simplify considerably. We find: 
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2 



(4.1) 



with the total conductance given by pi[| 



5Diff(e) 



2 



(4.2) 





tan(^) 


2 


'Ns 


cos(0(e)) 


2 




l+sin(e(£)) 






A(£)2 - 




2 



As discussed in the previous section, these two contri- 
butions to the Andreev conductance can be interpreted 
as representing two distinct physical consequences of the 
proximity effect. The direct term gDir represents pro- 
cesses in which the dot mimics a bulk superconductor, 
with Andreev reflections effectively occurring locally at 
the N-dot interface. It decreases monotonically with volt- 
age and magnetic field, going to zero at large voltages 
or magnetic fields. This reffects the fact that the in- 
duced superconductivity effect is sensitive to the aver- 
aged amplitude for Andreev reflection, which is largest 
dX V = B = 0. On the other hand, the diffusion term 
5Diff increases monotonically with V and _B, tending to 
the classical result for two conductances in series: 



g = 5Diff = 



(2jVg)(jVjv) 
2Ns + Nn ' 



(4.3) 



It represents a contribution from Andreev quasiparticles 
in the dot, and is thus sensitive to the dot's density of 
states. In the present case of ballistic contacts, the den- 
sity of states does not have any BCS-type peak and thus 
the diffusion term rises steadily with applied V ov B. 

Given that the direct and diffusion contributions react 
in opposite fashions to an increase in or i3, it is not sur- 
prising that a non-monotonic V or B dependence of the 
total conductance can be found if the relative strengths 
of these two terms are varied. The latter can be achieved 
by tuning the ratio of Nn/Ns- If Nn < Ns, the di- 
rect term will be dominant at V, i? = 0, and we expect a 
monotonic decrease of G as a magnetic field or finite bias 
are applied. In the opposite hmit ^ Ns, the diffu- 
sion term dominates at V,B — 0, and G is expected to 
increase with V or B. A non-monotonic V or B depen- 
dence can thus be anticipated in the intermediate regime 
where Nn and Ns are comparable. 

To quantify the competition between the direct and 
diffusion contributions, we examine these terms at ^ = 
B = 0. This is done by solving the self-energy equa- 
tions ( 2.15a| )-(2.17c) to determine 9{e), and then substi- 
tute this into Eqs. (4.1) and (4.2). Letting A^sum — 



^N% + QNnNs 



g-Dh = (-/Vsum 
ffDiff — (A/'sum 



A^J, we obtain: 

Ns- 
Ns 



N, 



N 



Ns,, 



N 



2Ns 



N 



(4.4a) 



Nn Asum 



Nn + Ns\ , , 
5Dir + 5Diff - {Ns ^Nn)\\- ) (4.5) 




FIG. 1^. Plot of the direct and diffusion contributions to the 
1/ = 0, -B = conductance as a function of Nn/Ns for trans- 
parent point contacts (Tjv — Ts = 1). At A^jv — fA'jv, the 
total V = 0, B — conductance is the same as the classical 

result Qclasa- 




B/B^ 



FIG. § 

Qclass ) I Qcla 



Plot of the excess conductance {g(B,V) — 
as a function of B at 1/ = 0, showing non- 
monotonic behaviour. From top to bottom, curves are for 
Nn/Ns = 0.39,0.45,0.54,0.68. The inset shows the excess 
conductance at _B = as a function of voltage; from top 
to bottom, curves are for Nn/Ns = 0.49,0.51,0.57,0.6. In 
general, re-entrance effects are maximized at Nn/Ns — 0.54, 
where the V = 0, B = conductance is the same as that large 
V or B. 
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Plotted in Figure]^ as a function ofNpf/Ns are the zero 
field values of guir , gms and 5, normalized by the high 
field conductance given by the classical formula (4.3). As 
Nn /Ns is increased from zero, the total conductance ini- 
tially decreases with ^Dir as the local Andreev reflection 
effect is suppressed, while at larger values it starts to in- 
crease with guis as the density of states at Ep returns to 
its normal state value. At Nn ~ 5-^5, we find that the 
total conductance at zero field is the same as the classical 
result. Not surprisingly, we find that re-entrance effects 
in both magnetic field and voltage are maximized here 
(see Figure H). 

It is instructive to make a comparison at this point 
to the non-monotonic re-entrance behaviour of diffusive 
NS systems (i.e. a diffusive normal metal in good con- 
tact with a superconductor). In the diffusive case, the 
field and high field conductances are the same, being 
equal to the normal state conductance. The lack of any 
change due to superconductivity at V — B = is usu- 
ally explained as an exact cancellation of the conductance 
doubling effect of Andreev reflection by a suppression of 
the density of states at the Fermi energy. A conduc- 
tance enhancement is however found at finite voltage, 
roughly at eV ~ Et = hD/L^, where D is the diffu- 
sion constant and L is the length of the normal metal. 
Quasi-classical calculations find maximum conductance 
enhancements on the order of 10% of the normal state 
conductance pl . 

In the present system, we can tune the relative signif- 
icance of the Andreev reflection enhancement and den- 
sity of states suppression terms by varying N^/Ns- We 
find that non-monotonic effects are maximized when we 
adjust this ratio to mimic the diffusive system, by in- 
sisting that the field and high field conductances are 
the same; this occurs at Nn ~ \/2Ns- The conduc- 
tance maximum in voltage is smaller however, being 
on the order of Q.QlGciassicaii and occurs roughly at 
eV — Es — Ns6/2tt. This is the inverse of the time 
needed to escape to the superconductor, and is thus the 
analog of Et in the diffusive system, which represents 
the inverse of the time needed for a particle to diffuse 
across the normal metal and reach the superconductor. 
Note also that if the normal lead was removed, the size 
of the induced gap A in the dot density of states is ^ Es- 

Finally, we also find pronounced non-monotonic be- 
haviour in magnetic field for this range of Nj^/Ns, with 
the magnitude of the effect being larger than that in volt- 
age (Fig. H). The maximum conductance occurs roughly 
at a flux $0" given by Eq. (3.7), the same flux that would 
be required to close the gap in the density of states in the 
absence of the normal lead. 



B. Tunnel Regime 

We turn now to the case where both point contacts 
contain opaque tunnel barriers (Ts,Tn ^1). In this 
regime, it is sufficient to consider the conductance to 
lowest non- vanishing order in T/v and Tg. Of particular 
interest here is the well-known "reflectionless tunneling" 
effect |2^]- the V = 0, B = Andreev conductance of a 
N-I-N-I-S structure (where / is an insulating region hav- 
ing transmission T) is found to be proportional to T, not 
to as one has for a single barrier. It is as though the 
Andreev reflected hole is not reflected at all by the tunnel 
barriers. A similar effect occurs for a N-I-S system where 
the normal region is sufficiently disordered. The origin of 
this striking behaviour is now understood to result from 
the constructive interference of trajectories which reflect 
many times between the two barriers, leading to a frac- 
tion ^ T of the conductance channels being open (i.e. 
having a transmission probability close to unity) p2| . 

In the present case, we are able to examine the effects 
of finite voltage and magnetic field on reflectionless tun- 
neling when a quantum dot separates the barriers. At 
large voltages or magnetic fields, the pairing angle 9{e) 
tends to its normal-state value of 0, and we find that 
the conductance is given to leading order by the classical 
series addition formula: 



5cla 



^NsTl + NnTn 



(4.6) 



For iVsTj < NnTn, this simplifies to c/ciass = ^NgT^. 
Unlike the case at V ~ 0,B = 0, there is no order T re- 
flectionless tunneling contribution here, as the necessary 
constructive interference is lost when electron-hole de- 
generacy or time-reversal symmetry is broken. Below we 
show that reflectionless tunneling does survive at small 
values of B and V, and describe how this contribution 
evolves as B and V are increased. 

We begin by solving Eqs. ( |2.17a| )-( |2l7c| ) for the self 
energy to lowest order in Tjv, Ts for B = 0; we find the 
pairing angle is given simply by: 



(e) = arctan 



Es 



E 



N 



(4.7) 



Es and E^ are the inverse escape times to the super- 
conductor and normal metal leads respectively, and are 
defined by: 



NnTn5 ^ NsTsS 



(4.8) 



Using this result for 9{e), we next write Eqs. (2.23a)- 
( p.231: ) for the conductance to lowest order in Ts,Tn- 
The direct contribution gY)ir corresponds to Andreev re- 
flection at the A^-dot interface and is order T^. The only 
order T contribution is found in the diffusion term goiff, 
which yields: 
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5'Diff(e) = NnTj 



EnEI 



N- 



E^e) y E^{e)+E% + El-e'^ 



where we have defined: 



o(t2), 



(4.9) 



(4.10) 



Eq. ( [4.9D gives the complete voltage dependence of the 
reflectionless tunneUng effect. At V — 0, it reduces to: 



E{e) = {El -E%-e^y+4 {EsEnY 



3Difi(0) 



{NnTn^ {NsTsf 



(4.11) 



which is similar to the formula found in |Q, general- 
ized to the case where the point contacts have different 
widths. Depending on the relative magnitudes of Eg and 
En, the conductance drops monotonically with voltage, 
or shows a maximum around eV ~ Es (see Fig. |^). 
The location of the maximum is at eV = \fTJ%Es if 
Eg 3> -Ejv- Similar behaviour is found for the magnetic 
field dependence of the conductance (sec inset of Fig. 0). 




FIG. ^. Conductance vs. voltage in the tunnel regime, 
Tat = Ts = 0.01, exfiibiting the enhancement of reflectionless 
tunneling at finite voltage. The parameter q is defined as the 
ratio Nn/Ns- The inset shows the V — conductance as a 
function of magnetic field for the same parameter choices. 

It is easy to understand the origin of this non- 
monotonic behaviour within our theory. When Ns be- 
comes larger than A^jv, the effect of the superconductor 
on the dot density of states becomes significant. As we 
have seen in Section 11, the induced density of states will 
have a sharp peak at Es when Is ^ 1; it is this peak 
that is manifesting itself in the conductance. 

It is interesting to note that Eq. (4.9) has the same 
form as what was found for a diffusive N-I-N-I-S system 
using a quasi-classical Green's function approach | |3^ . In 
that system, the quantum dot is replaced by a diffusive 



normal wire of length d, width W and mean free path 
I, and the normal and superconducting leads are not at- 
tached via point contacts, but are also wires of the same 
width. In the limit T/v,Ts <C l/d the resistances of the 
tunnel barriers dominate, and an expression identical to 
Eq. (4.9) was obtained, but now the energies En, Es are 



given by: 



E 



N 



Tivp 
Tn^.Es 



Ts 



(4.12) 



where is the Fermi velocity in the normal inter-barrier 
region. This is identical to the definition of E^ and Es 
in Eq. (4^) if one takes S to be the level spacing in 
the wire. Note that this correspondence is not surpris- 
ing; previous studies have also found that adding strong 
tunnel barriers makes a diffusive normal wire with many 
channels equivalent to a quantum dot. 



V. CONCLUSIONS 

We have studied the voltage {V) dependence of the 
Andreev conductance of chaotic and integrable quantum 
dots; we have also examined the magnetic field {B) de- 
pendence in the chaotic case. In the regime where the 
contact to the superconductor dominates, we find that 
the voltage dependence of the Andreev conductance is 
extremely sensitive to the nature of the dot's classical 
dynamics- in the chaotic case, the dot itself mimics a 
superconductor and the conductance is initially flat in 
voltage, whereas in the integrable case the conductance 
falls off linearly with voltage. This effect is large in that 
the entire conductance forms the signature of chaos vs. 
integrability; also, it does not require any ensemble aver- 
aging to be done. Both these facts make it particularly 
amenable to experiment. Also important in this regard 
is that one does not require an extremely clean contact 
between the dot and the superconductor as long as the 
contact is large. 

We have also studied non-monotonic "re-entrance" 
phenomenon in the V and B dependence of the Andreev 
conductance of a chaotic dot. We find that such be- 
haviour is ubiquitous, and is the result of two competing 
processes- Andreev reflection at the dot-normal lead in- 
terface, which decreases asV or B increases, versus quasi- 
particles being injected into the dot, which increases with 
V and B. 
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APPENDIX A: CALCULATING THE ANDREEV 
CONDUCTANCE FOR B 7^ 



In this appendix, wc outline the method used for cal- 
culations at non-zero magnetic field B. Formally, the 
magnetic field enters the model in a different fashion 
than a voltage difference. The latter is dealt with by 
the fact that the Andreev scattering matrix S defined in 
Eq. (2.12) is an explicit function of energy. The field 
dependence however does not appear directly in the ex- 
pression for S, but rather only emerges in the averaging 
procedure- as we use the Pandey-Mehta distribution de- 
fined in Eq. (2_^), the ensemble of random matrices is 
itself a function of field. 

As in the calculation at B = 0, the first step in obtain- 
ing the conductance is to calculate the averaged matrix 
Green function (Q) defined in Eq. (2.14a). The diagrams 
used arc the same as the B — case, but now use of the 
distribution ( ^.2|) leads to a different self-energy S in the 
Dyson equation: 



J 



y2 



A2 



- 1 = - 



1 

M 



S(£,7) = 1m ® 



1m (S> 



See (1 - 2j)T,eh 

(1 - 27)S,ie 

f (tr Gee) 

M W27-l)(tr Gh 



(Al) 



(27-l)(tr Geh) 
(tr Ghh) 



where 7 is a function of magnetic field (see Eq. 2.3). The 
additional factors of (1 — 27) here reflect the fact that 
breaking time-reversal symmetry suppresses off-diagonal 
superconducting correlations. Solving the Dyson equa- 
tion ( ^.14a ) to leading order in 1/M, we find that rela- 
tions ( |2.15a ) and ( ^.15b| ) continue to hold, meaning that 
we may still parameterize the self-energies in terms of a 
pairing angle 9{e,^) through Eq. ( p. 16a ). The equation 
determining ^(^,7) now takes the form: 



tan(0(e,7)) 



7V5Q5- 27 sin(0) 



N 



(A2) 



'2<5 



where Qn a nd Qs are functions of 9 defined in Eqs. 
(^.17b| ) and (|2.17c| ). The l/M corrections to the self- 
energy in the presence of a magnetic field read: 



ITS 



NsQs (sin(^^) + Zs) + NnQn (cos(0) + Zn) + + 2-fsm'{9) 

Zd 



(A3) 



r 



The next step in the calculation is to sum diagrams for (^.23b| ), as factors of (1 — 27) now appear in graphs where 

the conductance. The necessary diagrams are the same the particle-hole indices of upper and lower branches do 

as those retained in the B = calculation (i.e. direct not match. These factors lead to 1/M corrections both 

and diffusion terms), although their evaluati on is differ- to the matrix inverse that arises when summing the dif- 

ent. The direct term is still given by Eq. ( 2.23a ), but fusion ladder, and to the matrix prefactors of the ladder, 

with 9{s, 7) now determined from Eq. ( A2). As discussed earlier, such corrections are important when 

The diffusion term acquires a form different from Eq. calculating the conductance. The result is: 
I 



N 



Q 



(1 - Ziyk\ sin(g)p + (Hat + 70^)1 sin(g)p + Hgl cos(g)|- 



where 
^b{9) 

Yn{9) 



2(^(1 + Zn^ (1 + 2Re[sin2(6l)]) + 4Zjv(l + Zjv)Re[cos(6')] - 
A{9) - 127sin2(6'), 

ni9) + 47 (1 - \cos{9)f + 4Re[sin2(6l)]) A(6l,7) + 87 (^Nn 
Zn 



{\cos{9)\^ - 2Re[sin2(6i)]) ), 



Z 



N 



Yn{9)+Ns 



Zs 



Ys{9) + -iY2{9) 



- 1 



1 - |cos(6l)|^ +4Re[sin^(6l)]) - 2ZAr (5Re[sin2(6l)]Re[cos(6')] + Im[sin2(6')]lm[cos(6')] 



+Z%{l - \cos{9)\^ - 2Re[sin2(6l)] (n- 2 |cos(6')p 

- |1 - Zssin(6i)|^^ (^1 - |cos(6')|V4Re[sin2(e')]) - Z||sin(6l) cos(6l)|^ 
= 1 - |cos(6')p (1 + 4Re {sin{9f)) + Re [3sin2(6') + 4sin'^(6')] + 4 |sin(6')|* . 




(A4) 

(A5) 
(A6) 

, (A7) 



(A8) 
(A9) 
(AlO) 
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Here A {0), n je), Un{9) and lls{0) are given by Eqs. 
( 2.22aD - ( 2.22d ). The above equations, together with Eq. 
( A2) for 0{e,-f), determine the Andreev conductance for 
arbitrary voltage and magnetic field. 
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